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Symmetries represent a fundamental constraint for physical systems and relevant new phenomena 
often emerge as a consequence of their breaking. An important example is provided by space- 
and time-translational invariance in statistical systems, which hold at a coarse-grained scale in 
equilibrium and are broken by spatial and temporal boundaries, the former being implemented 
by surfaces — unavoidable in real samples — the latter by some initial condition for the dynamics 
which causes a non-equilibrium evolution. While the separate effects of these two boundaries are well 
understood, we demonstrate here that additional, unexpected features arise upon approaching the 
effective edge formed by their intersection. For this purpose, we focus on the classical semi-infinite 
Ising model with spin-flip dynamics evolving out of equilibrium at its critical point. Considering both 
subcritical and critical values of the coupling among surface spins, we present numerical evidence of 
a scaling regime with universal features which emerges upon approaching the spatio-temporal edge 
and we rationalise these findings within a field-theoretical approach. 

PACS numbers: 64.60.De, 64.60.Ht, 68.35. Rh 



Introduction - Thanks to the advances in miniaturisation of the past decades, de- 
vices have reached a size at which boundary effects cannot be neglected; in parallel, 
the physics of surfaces and interfaces has attracted an increasing interest, in particular 
concerning non-equilibrium dynamical processes, as many applications involve sudden 
changes in control parameters. Typically, describing these processes requires knowledge 
of the many microscopic details that vary widely from system to system. However, in suit- 
able circumstances, collective phenomena emerge which make only few coarse-grained, 
mesoscopic properties relevant for the dynamics: in fact, it is now well-established, both 
theoretically and experimentally, that the behaviour of a statistical system close to a 
continuous phase transition can be characterised by few quantities, such as exponents 
and scaling functions, which depend only on the range and symmetries of the underly- 
ing interaction and on the dimensionality of the space. All the microscopically different 
systems sharing these same gross features form the so-called universality class of the tran- 
sition. Within each of these classes, the various thermodynamic and structural properties 
typically show, in the neighbourhood of the critical point, leading algebraic behaviours 
characterised by common exponents, which constitute the hallmark of the transition. 
In turn, upon approaching it, the relevant contribution to the various thermodynamic 
quantities is effectively determined by the fluctuations of the so-called order parameter 
ip (e.g., the local magnetisation for an Ising ferromagnet). 
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The emergence of universality is currently understood within the framework of the 
renormalisation group (RG) [l|, 0|, which is among the most important theoretical 
achievements of the past forty years. RG transformations highlight scale invariance, a 
feature that emerges upon approaching a continuous phase transition, as the only meso- 
scopic length-scale present — i.e., the correlation length £ of the fluctuations of the order 
parameter — ideally diverges; because of it, one can focus on effective, coarse-grained 
descriptions and models which retain only the gross features of the original system and 
which are effective at space and time scales larger than the microscopic ones. Originally 
devised for describing the local behaviour of finite samples far from the boundaries (i.e., 
in their bulk), this theoretical tool has been subsequently extended in order to account 
for the presence of surfaces (e.g., see Refs. 0,01)- Their introduction alters locally the 
properties of the system and breaks translational invariance: within the RG approach, 
a new set of relevant parameters has to be introduced in order to account for the gross 
features of the boundary. As a result, novel singularities might emerge upon approaching 
the surface, which split the original universality class in surface subclasses characterised 
by a set of boundary exponents and scaling functions associated, e.g., with the algebraic 
behaviour of the cp correlation functions in its proximity 0, 0| . In general, these expo- 
nents cannot be inferred from the bulk ones. A number of analytical [5|-|7j], numerical 
0, 0| and experimental (see, e.g., Ref. [Io| ) studies investigated primarily semi-infinite 



and film geometries, whereas wedges, edges 111 |12|. as well as curved and irregular sur- 
faces 0,[l3| were studied to a lesser extent. Universal features emerge also in the dynamic 
behaviour at equilibrium (both in infinite and finite systems) 14j . and the universality 
class is further split depending on the gross features of the dynamics, such as the possible 
global conservation of the order parameter fl5. 16|. 



A sudden quench of the system's temperature breaks time-translational invariance — 
a characteristic symmetry of stationary states — and, analogously to the case of spatial 
surfaces, effectively introduces a sharp temporal boundary for the evolution and causes 
non-equilibrium behaviours to emerge [l7, 18 1. In particular, the memory of the initial 



state naturally affects the dynamics close to this boundary and, in certain circumstances, 
it is also responsible for the later occurrence of ageing phenomena (rfl - tlij . Within this 
picture, the dynamics at "short" times becomes akin to the behaviour of static quantities 
at short distances from a spatial surface [l6| and therefore RG techniques can be applied 
in order to study the short-time and ageing behaviours and their emerging universal 
features [3, [lj|. In order to investigate the subtle interplay between the breaking of 
space- and time-translational invariance we study here the semi-infinite Ising ferromag- 
net quenched from a disordered state to its critical temperature at time t = 0. This 
relatively simple, but paradigmatic case can be used to verify whether effects beyond 
those resulting from each separate breaking [2(| emerge. We use both numerical and 
analytical methods to show that, unexpectedly but similarly to the case of a spatial 



wedge [ll|], a so- far undetected power-law behaviour described by a critical exponent 9e 
emerges upon approaching the effective edge formed by the intersection of the spatial 
and temporal boundaries. 
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Scaling behaviour - Analogously to the Ginzburg-Landau (GL) theory of phase tran- 
sitions, the collective dynamical properties of Ising ferromagnets can be studied in terms 
of a continuum field theory for the order parameter ip(x,t) (which represents here the 
coarse-grained local magnetisation) , the thermal stochastic dynamics of which is captured 
by a Langevin equation (l5| 

where rj(x,t) is a white Gaussian noise with zero mean and variance (rj(x,t)r](y, s)) = 
2J7fcsT(5 (x — y) 8 (t — s) and H is the static GL free-energy Q| 



H[(p]= J d d x 

{i±>0) 



(Vy>) 2 r 2 g 4 c 2 
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(2) 



Here x± stands for the coordinate normal to the surface located at x± = 0, d is the 
space dimensionality and, denoting by T the temperature and by T c its critical value, 
r oc T — T c . Whereas g > characterises the interaction strength in the bulk, the so- 
called surface enhancement cq is related to the difference between the values ,J\, and J s of 
the coupling constant among the spins in the bulk and within the surface, respectively. 
Depending on its value, one distinguishes three cases 0, Q|: (i) cq —> +oo (ordinary 
transition) corresponds to J s <C Jt which implies that the surface ordering is induced by 
the bulk one; (ii) cq — > — oo (extraordinary transition) corresponds to the case J s 3> Jb, in 
which the surface acquires a non- vanishing magnetisation independently of the bulk; these 
two instances are separated by a critical value (iii) cq = c sp (special transition), which 
corresponds to a certain ratio J s /Jb (—1.5 for the three-dimensional Ising model (2l|). 
The last two cases are present only for d > 3. Hereafter (■) indicates the average over the 
possible realisations of the thermal noise 77 and of the initial state; for convenience, we set 
the scales of time and T such that n = fc^Tc = 1 and fix the system at its bulk critical 
point T = T c . The initial condition is chosen to be an effective high-temperature state 
both in the bulk and at the surface, characterised by local correlations (tp(x, t = 0)ip(y, t = 
0)) = Tq 5(x — y) with tq much larger than any other scale, i.e., we will actually consider 
the limit To — > +00. In the following we shall concentrate on the ordinary and special 
transitions. The corresponding average order parameter {<p(t)) was studied in Ref. (2o| . 
where a power-counting argument was proposed in order to argue that no new field- 
theoretical divergences arise at the spatio-temporal edge (t = 0, x± = 0). Accordingly, 
the exponents of the associated algebraic singularities should be a combination of those 
characterising the static bulk (/3, v) and surface {fix) behaviours and the bulk equilibrium 
(z) and non-equilibrium (9) dynamics. For the total magnetisation <&(xx,t) of the plane 
at a distance x± from the surface (i.e., the spatial integral of ip(x,t) along the d — 1 
directions parallel to it) the two-point and two-time correlation function 



C(x,t;y,s) = ($(M)$(y,s)) 



(3) 
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is expected to scale as 



xFc ((^r, ( 4 ) 

\ t-s £ — s £/ 

where A is a non-universal constant which makes the scaling variables dimensionless and 
is of order unity for a suitable choice of time and space units. Fc is a scaling function 
(depending, inter alia, on the surface transition) with finite Fc(0, 0, 0) 7^ and such that 
the usual equilibrium scaling is recovered in the bulk and at the surface when t — >• s. 



For the d = 3 Ising class, a = (d - 1 - 2f3/v)Jz ~ 0.96 [23j, |2J, 1/ ~ 0.63 [23J is the 



exponent of the correlation length, j3 ~ 0.33 [23| and /3i l0 rd — 0.80 (or /3i )Sp ~ 0.24) 



13j are the order parameter exponents in the bulk and at the surface, z ~ 2.04 [24 1 



is the bulk dynamic critical exponent, whereas 6 ~ 0.15 [24| governs the short-time 
behaviour in the bulk. Note that in eq. Q /?i is the only exponent that distinguishes 
the ordinary from the special transition. We now define what we refer to as the "edge 
regime" (E): x, y and t are fixed such that f « ( - s, while s varies within the range 
f < s < t; correspondingly, eq. Q yields C(. . . , s) oc s 1-61 , where the proportionality 
constant depends, inter alia, on ( > s, x z /t and y 2 /i <C 1. Note that, according to 
eq. ([J), the assumption y z <C s is not necessary for observing this power law; however, 
it becomes important for the extension discussed further below. Conversely, we name 
"short-time" (ST) the regime in which s is much smaller than any other scale (which 
can be obtained from the edge regime by moving s to the domain s <C y z )- In summary, 
eq. @ predicts that the power-law behaviour of C as a function of s in the edge regime 
is the same both for the special and the ordinary case. 

Numerical analysis - Using the standard setup for studying surface criticality [l3j|, 
we simulate the Ising model on a three-dimensional cubic lattice made up of H = 40 
consecutive planes with 60 x 60 spins per plane. Within each of these planes, periodic 
conditions are imposed at the boundaries. The coupling constant between two surface 
spins is J s , while in the rest of the lattice it is set to J& = 1. The ordinary case is 
realised for < J s < 1.5, whereas the special transition is known to occur at J s ~ 1.5 



2l| . The system is prepared at t = in a completely disordered state corresponding 
to infinite temperature T and it subsequently evolves at its bulk critical value T c = 
4.5115J(,/fcs with Glauber dynamics. One time step corresponds to one sweep where 
on average every spin of the lattice has been updated once; as we are interested in the 
behaviour close to the temporal surface, only rather short simulation times are needed 
and therefore we can restrict ourselves to rather small system sizes. In order to rule 
out the influence of finite-size effects we also analysed some larger lattices. For the 
purpose of investigating the edge regime we study the two-point correlation function 
C(x,t;y, s) defined in eq. ([3]); here, x = 1 corresponds to the surface, whereas x = H/2 
to the midplane of the lattice. The statistical average is taken over a large number of 
independent runs with different realisations of the thermal noise and of the initial state. 
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FIG. 1: (Color online) Time dependence of the correlation function C(x,t;y,s) of the plane 
magnetisation of a three-dimensional Ising model evolving with Glauber dynamics after a quench 
from a completely disordered state to the bulk critical temperature. The bulk-surface correlation 
function C with x = H/2 = 20 and y = 1 is plotted as a function of s for t = 10 3 sweeps 
(corresponding to the "edge regime" discussed in the main text). Data points result from 
averaging over 5 x 10 s and 3 x 10 6 independent runs for the special (red upper dots) and ordinary 
(blue lower dots) case, respectively. The black solid lines correspond to power laws with the 
exponents reported to the right; the parallel dashed lines, instead, indicate, for comparison, the 
power law with exponent 0.85 observed within the short-time regime (see fig. 0) and predicted 
by eq. ©. The data points with s < 6 have been disregarded as they might be affected by 
microscopic effects. 



In fig. Q] we show the bulk-surface correlator C(x = H/2,t;y = l,s) for the ordinary 
case J s = 1 and the special one J s = 1.5, where we fix t = 10 3 and vary s within a 
range that, under the assumption A ~ 1, corresponds to the edge regime (we verified 
that different values of t give the same result). We observe that C(. . . , s) ~ s p increases 
algebraically, but with an exponent p sp = 0.71(2) for the special case and p or d = 1.05(2) 
for the ordinary one. This unexpected result p sp ^ p or d clearly shows that different 
surface universality classes display a different behaviour at the spatio-temporal edge, 



while eq. ([4]) requires p sp = p orc i = 1 — 9 ~ 0.85 [24 1. 



Field-theoretical results - In order to rationalise these findings, we carried out a field- 
theoretical perturbative calculation of C up to the first order in e = 4 — d (d = 3 in 
the simulations) ; beyond correctly reproducing the previously-known results concerning 
separately each of the two boundaries x± = 0, t — 0, this approach highlights indeed 
the emergence of an additional dimensional pole cx e _1 Q when the coordinates of the 
correlation function §3§ are fixed at the spatio-temporal edge y = s = 0, analogously 



to what happens for the static critical behaviour in a wedge [ll| . This pole is actually 
related to a dangerously irrelevant edge operator ipg, which appears as the most relevant 
field in the short-distance expansion Q of the order parameter <p close to the edge; 
in the present case, due to the properties of the boundaries, ifE can be identified with 
dt<p(t. x±) \ t _ x _ in the special case and with dtd x± ip(t, x±) \ t _ _ in the ordinary one 
[i, 18 1 . In order to account for its presence, one has to introduce a "radial" coordinate 
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FIG. 2: (Color online) Time dependence of the correlation function C(x = 20, t = 10 3 ;y,s) for 
the special and ordinary transitions and various values of the position y. For the 5 lowermost 
(uppermost) curves, corresponding to the ordinary (special) transition, y increases from bottom 
to top (top to bottom) and takes the sequence of values reported for the ordinary case to the 
right. This plot highlights the crossover from the "edge" behaviour (thick solid lines), observed 
for i/ z < s in accordance with fig. [T] to the "short-time" behaviour for y z ^> s, indicated by 
the thick dashed line, which is common to the two instances. For the ordinary case and y — 4 
this crossover can be observed as a function of s, as highlighted by the thin solid and dashed 
lines that meet at s ~ (Ay) z , which yields the indicative estimate A ~ 2. The details of the 
simulations are the same as for fig. [T] 



(Ay) z + s, with A fixed as in eq. which controls the distance from the edge. Because 
of ifE, the scaling form (U) acquires a new algebraic term, which can be made explicit 
by the substitution 

where now Fc(0, 0, 0) ^ is finite. While, as expected, this new scaling form reproduces 
the results of eq. Q when approaching either the spatial (e.g., fixed s and y — s- 0) or 
the temporal surface (fixed y and s — > 0), the additional term in eq. ((5]) affects the 
behaviour of the correlation function in the edge regime as C(. ..,s) ~ s 1 " 6 "^, which 
yields p = 1 — 8 — Be- The critical exponent Be encodes the anomalous dimension of ipE 
and, up to O(e), turns out to be jiil 



'E,sp 



V3-l)± and ^- = -(l-^ W 



for the special and ordinary case, respectively. Table U presents the comparison between 
the numerical estimates of p from fig. Q] and the corresponding analytic expressions re- 



sulting from Eqs. (0} and ©, with (up to 0(e)) B = e/12 [18fl and B E as in eq. © 



specialised to the three-dimensional case e = 1. Even if lacking quantitative accuracy, 
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Table I: Comparison between the numerical and analytical estimates of the exponent p which 
controls C(. . . , s) ~ s p in the edge (E) and the short-time (ST) regimes according to the scaling 
forms Eqs. and (J5]). 





Monte Carlo 


Analytical (e 


= 1) 




Figs. [1] and EH 


eq. g]) 


eq. © 


E, sp 


0.71(2) 


0.917 


0.856 


E, ord 


1.05(2) 


0.917 


0.952 


ST, sp & ord 


0.85(2) 


0.917 


0.917 



the additional term introduced in eq. (|5|) correctly captures the qualitative behaviour as 
a function of s of the correlation function C(x, t\ y, s) near the spatio-temporal edge: in 
the case of the ordinary and special transition C(. . . , s) grows respectively faster and 
slower than it does in the short-time regime. Figure [2] presents a study of the crossover 
between the edge regime of fig. Q] and the short-time regime, which occurs upon increas- 
ing y above the scale set by s 1 ^ /A. This crossover is properly captured by the scaling 
function in eq. ([5]) because the additional multiplicative factor becomes independent of 
s for s <§; t, so that C(. . . , s) ~ s 1-6 \ as predicted also by eq. (j4|. Conversely, for a fixed 
y =/= 1, the crossover between the short-time and the edge regime occurs upon increasing 
s above (Ay) z . With the present data, we could find clear evidence of this crossover for 
the ordinary case with y = 4 (see fig. [2]). Note that the edge regime is not the only 
one affected by the aforementioned factor; however, we focused on it in order to test a 
qualitative difference between the predictions of eq. (|4]) and ([5]) , which does not hinge on 
a difficult, quantitatively accurate analytic determination of 0e- 

Conclusions - By studying the non-equilibrium dynamics of the Ising model in prox- 
imity to a surface, we have identified a regime in which the spatio-temporal edge, defined 
as the intersection between the spatial and temporal boundaries, affects the scaling be- 
haviour of the correlation function; we rationalised and qualitatively accounted for the 
numerical evidence within a perturbative, field-theoretical framework. The consistency 
of these approaches, which concern two completely different systems at the microscopic 
level, i.e., a discrete spin model on a lattice and an interacting field theory in the con- 
tinuum, provides additional support to the universality of the edge behaviour, which 
could be expected as a consequence of scale invariance at the critical point. This ac- 
tually indicates that the scaling near the edge investigated here depends only on the 
gross features (symmetry, dimensionality and short-range interaction) which are shared 
by both models. Analogous results have been obtained for the linear response function 
and within the vector generalisation of the present model It would be desirable 

to extend the present investigation to different static, dynamic and surface universality 
classes of experimental relevance. While the first experiments probing the static surface 
and equilibrium dynamical bulk properties of materials date back to the 1970s 
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techniques with the sufficient accuracy to study bulk non-equilibrium behaviours have 
been available only since the 1990s [27j]. Surface dynamics in condensed matter systems, 
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instead, have not been observed until recently, though with a different purpose and not 
in the critical regime [2^. Extending recent investigations of ageing phenomena in liq- 
uid crystals (2^ to the proximity of surfaces might provide a viable alternative for the 
experimental test of the present predictions in systems undergoing an Ising transition. 
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